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Abstract 

In this paper, we study some important statistics of the random graph in the Buckley-Osthus 
model, where t is the number of nodes, kt is the number of edges (so that k E N), and a > is the 
so-called initial attractiveness of a node. This model is a modification of the well-known Bollobas- 
Riordan model. First, we find a new asymptotic formula for the expectation of the number R(d,t) 
of nodes of a given degree d in a graph in this model. Such a formula is known for a 6 N and 
d £i/ 1 0°O+ 1 ). Both restrictions are unsatisfactory from theoretical and practical points of view. 
We completely remove them. Then we calculate the covariances between any two quantities R(di,t), 
R(d2,t), and using the second moment method we show that R(d,t) is tightly concentrated around 
its mean for every possible values of d and t. Furthermore, we study a more complicated statistic of 
the web graph: X(di,d2,t) is the total number of edges between nodes whose degrees are equal to d\ 
and c?2 respectively. We also find an asymptotic formula for the expectation of X{d\,d2,t) and prove 
a tight concentration result. Again, we do not impose any substantial restrictions on the values of 
d\,d2, and t. 



1 Introduction 



> 

in 
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The real world has many interesting structures which can be thought of as graphs. A typical example 
^ is the World Wide Web: one can consider web pages to be nodes of a graph and hyperlinks to be edges. 

i One of productive methods for studying these graphs involves investigation of a suitable random graph 
j>! hiodel. 

First models of random graphs were constructed and investigated long ago. Classical models and 
results are systematized, for example, in [I] and [IT] . However, they are not suitable for approximation 
. 5^ of dynamically changing and non-uniform networks. In particular, the degree sequences of the graphs in 
these models are very far from those observed in reality. 

Recently other models of random graphs were constructed to more closely match the growth of real 
networks. One of the first descriptions of such a model belongs to the article [3] by Barabasi and Albert. 
The authors of this article introduced the "preferential attachment" rule. Models following this rule assign 
the probability of a new edge to a node according to the current degree of this node, so more "popular" 
nodes are more attractive for new edges. 

However, the article [3] did not contain a precise model, leaving some parameters unspecified. Variations 
of these parameters can significantly change properties of arising graphs, as shown in [5], so one needs 
something more explicit for theoretical investigations. Bollobas, Riordan et al. proposed an explicit model 
in [6] based on the preferential attachment rule. In the same article, they rigorously proved a theorem 
concerning the degree sequence of a graph in this model. Namely, they showed that the number of nodes 
with degree d in their model decreases proportional to d~ 3 . The same quantity in real networks decreases 
proportional to <i~ 7 with different 7 for different networks, following the so called "power law" . 
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The Bollobas-Riordan model has only one parameter, a natural number representing the ratio of the 
number of edges to the number of nodes. Thus, on the one hand, the Bollobas-Riordan model does 
certainly match some real networks by explaining the power law. But, on the other hand, the number of 
parameters in this model is small and does not allow to obtain the power law with an exponent, which is 
not equal to —3. 

In the Bollobas-Riordan model, the probability that a node is target for a new edge is proportional 
to the degree of this node. In [S] and jH] two groups of researchers independently proposed to add to 
the model one more parameter — an "initial attractiveness" of a node which is a positive constant not 
depending on the degree. Equivalently, the probability in the proposed model is a linear function in the 
degree. However, in the papers [8] and [9], we find only some heuristic arguments. 

In [7] Buckley and Osthus gave an explicit construction of the above- described model and rigorously 
proved a theorem concerning the degree sequence of a graph in this model when all the parameters are 
natural numbers. 

Among many articles in this area, we also quote |12j . The model investigated in this article differs from 
the Buckley-Osthus model, but the difference is small, so the results are comparable. The article deals 
with the case when parameters are not necessarily natural. However, the proven theorem only works for 
fixed degree d when the number of nodes tends to infinity; Bollobas et al. as well as Buckley and Osthus 
allowed d to grow as some small power of the number of nodes. 

There are many other random graph models intended to approximate real networks. We refer the 
reader to [5] and [10] for surveys of such models and corresponding results. 

We study the Buckley-Osthus model of a random graph. Our first goal is to give a significant im- 
provement of the above-mentioned theorem from [7J using a completely different method. We find an 
asymptotic formula for the expectation of the number of nodes with degree d without any upper bound 
on d and with an estimation of the error term. We also prove a tight concentration result. 

Since the Bollobas-Riordan model is a special case of the Buckley-Osthus model, our results are also 
applicable to it. So, again, we get a substantial improvement of the main theorem from [6]. 

Our second goal is to study the following quantity. We fix two numbers d\ and d 2 - We consider a node 
with degree di and a node with degree d 2 - Then, we calculate the number of edges between these nodes. 
When there are several choices for nodes of given degrees, we calculate the mean value. Since the number 
of nodes with a fixed degree is known to have tight concentration around its expectation, it is sufficient to 
examine the total number of edges linking a node with degree d\ and a node with degree d 2 - Here we also 
obtain an asymptotic formula for the expectation and prove a tight concentration result. 



2 The model and formulation of results 

The Buckley-Osthus model has two parameters, a natural number k and a positive real number a. 
The number k is the ratio of the number of edges to the number of nodes. We assume that a and k are 
constants, so by default all other constants may depend on them. The Bollobas-Riordan model is a special 
case of this model with a = 1. 

The model is defined in two stages. At the first stage, a probability space is constructed. The 

elements of are undirected graphs with nodes represented by numbers 1, . . . , t and with t edges. The 
space H^l contains only one graph with one node and one loop. Any graph in is obtained from a 
graph in by adding a new node t and a new edge between t and a node 7 G {1, . . . , t} so that 



Pr(7 = s) 



l(sM+a l<s<t-h 



(a+l)i-l ' x ^ ° ^ " x > 
S = t, 



(a+l)t-l' 

where deg t-1 denotes the degree of a node in the graph from H^ x \ At the second stage, a final probability 

space H® k is constructed from as follows. We take any graph from H^' . It has kt nodes and kt 

edges. We identify the nodes 1, . . . , k; k + 1, . . . , 2k; . . . obtaining t new nodes, and we keep all the edges 
obtaining multiple edges and even multiple loops. 



We study the number of nodes of degree d in \ as a function of d and t. We denote this random 
quantity by R(d,t) and the value of its expectation by r(d,t) = ER(d,t). 

If d < k, then clearly R(d, t) = 0, so it suffices to study the case d ^ k. We start by considering r(d, t). 

Theorem 1. Let d^k. The expected value of R(d,t) is 

. , , B(d-k + ka,a + 2) „ (\ 

r * = fri — 1 + 0a . fc w 

B(/ca, a + 1) \« 

T/ie asymptotic behaviour of the coefficient when d grows is 

B(d-k + ka,a + 2) ^ T(a + 2) d _ 2 _ a = T(ka + a + 1) ^__ 2 _ Q 

B(ka,a + 1) ~B(ka,a + l) l ° ' T(ka) 

A similar result was obtained in [ij (with some factorials instead of Gamma- and Beta-functions). 
However, for that result, it was essential that a e N and d ^ ^i/ioo(o+i)_ Another result, which can be 
compared with the one of Theorem 1, is proved in [12] . It concerns a bit different model, but, nevertheless, 
it is rather close to our investigations. In this result, a can be any positive real (and so analogous Gamma- 
and Beta-functions appear in its statement). However, its proof essentially uses the assumption that d is 
just a constant. In our Theorem 1, we do not have any restrictions on d and a, and we use a completely 
different method to prove it. 

In fact, Theorem 1 gives an entire picture of what happens to the quantity r(d,t). If d = o fi^+H, 

then Theorem 1 yields the main term of r(d,t). If d = f2 (t^j, then r(d,t) tends to zero as t — > oo, 

which means that with high probability there are no nodes of degree d in a graph in the model. 
Now we want to study in detail the quantity R(d,t). 

Theorem 2. Let d\^k and di ^ k. The covariance between R(di,t) and R{d,2,t) is 

cov(R(d h t),R(d 2 ,t)) = O a)k {{d^ 2 - a + df~ a ) t + d^d^ 1 ) . 

Substituting d\ = d 2 = d in Theorem 2 and using Chebyshev's inequality, we obtain the following 
result. 

Corollary 1. If d — d(t) ^ k and ip(t) — > oo as t — > oo, then 

B(d — k + ka, a + 2 



< (yd- a -H + cT 1 ) i/j(t) (1) 



B{ka, a + 1) 

with probability tending to 1 as t —t oo. 

Let us discuss the meaning of Corollary 1. 

When d ~ Ct«+2 with some constant C, both r(d,t) and Vd~ a ~ 2 t + d^ 1 are 0(1). For smaller values 
of d (i.e., d = o U*+ 3 ^), inequality (JT]) implies the equivalence (with probability tending to 1 as d, t — > oo) 

(a + l)r(fca + « + l) 
R(d,t) ^ d t. 

For larger values of d (i.e., t°+2 = o(d)), inequality (P) means that R(d,t) = o(l). Since R(d,t) is an 
integer number by definition, R(d,t) = (again, with probability tending to 1 as d, t — > oo). Thus, we 
have an almost entire picture of what happens to R(d,t). 

We also study the total number of edges linking a node with degree d\ and a node with degree d 2 - We 
denote this random quantity by X(d 1 ,d 2 ,t). When d\ = d 2 , we count every edge twice, but do not count 
loops. 



Theorem 3. Let d\^ k and d 2 ^ k. There exists a function cx{di, d 2 ) such that 

EX(d 1 ,d 2 , t) = c x (d 1 , d 2 )t + O a>fc (l) 

and 

Y{di - k + ka)T(d 2 - k + ka)Y(di + d 2 - 2k + 2ka + 3) 

u 2 > = r(di - k + ka + 2)T(d 2 - k + ka + 2)T(d 1 + d 2 - 2k + 2ka + a + 2) X 

, T(ka + a + 1) / „, , ,. (di - k + ka + l)(d 2 - k + ka + 1) 

x ka(a+ 1) v l + 0(di,d 2 ) 



r(ifea) V J (d 1 + d 2 -2k + 2ka + l)(d 1 + d 2 -2k + 2ka + 2)J ' 

2 . * T(A;a + l)r(2A;a + a + 3) 

" 4 + TT^ ^ ^ ^ °r(2A;a + 2)r(A;a + a + 2)- 
When both d\ and d 2 grow, the asymptotic behaviour of cx is 

1 d\d 2 

c x (di, d 2 ) = ka(a + I)——— v - , 9 I 1 + U„ k I — + — + 



(A A\ u i i 1 ^ r(fca + a + l)(d l + d 2 ) l - a / /l 
: (d 1 ,d 2 ) = Jba(a + l)^^ 5 ^^—^l + O^ 



d 2 (di + d 2 ) 



Note that the last formula in Theorem 3 does not give an asymptotic behaviour if d\ and d 2 grow so 
that ^ tends to a finite nonzero limit. The precise bounds show that the term ^ dl+ ^p — still gives the 

correct order of growth for cx, but the coefficient can differ from ka(a + 1) ^^°^ . And in fact, the 
coefficient differs. 

Theorem 4. Let d\, d 2 ^ k, c > 0. Then 

P (\X(d 1: d 2 ,t) - EX(d u d 2 ,t)\ ^ c(di + d2)Vkt) <2exp^-^. 



In particular, if c(t) — >■ oo as t — > oo, then \X — EX\ < c{t){d\ + d 2 )ykt with probability tending to 1. 

From Theorems 3 and 4, we immediately obtain the following assertion. 
Corollary 2. // (d\ + d 2 ) a d\d\ = o (Vi), then with probability tending to 1 as d 1: d 2 , t — > oo 

X(di ) d 2) i) ~ c x {di,d 2 )t. 



The mean value of the number of edges between one node with degree d\ and another node with 
degree d 2 is R ^f^^l ^ • Since the quantities R(d,t) and X(d 1 ,d 2 ,t) are tightly concentrated around their 
expectations, the main term of the ratio is 



T(ka + 1) dfd§(di + d 2 ) 



1-a 



(a+ l)T(ka + a + 1) 

cr\ +Vta + 

di 



Again, the constant factor can differ if d\ and d 2 grow so that ^ tends to a finite nonzero limit, but the 
order is correct even in this case. 



3 Proof of Theorem 1 

For a property P, we denote 

[P] 



P holds, 
otherwise. 



First of all, we reformulate the model without references to H^\. The probability space H^l obviously 
consists of one graph with 1 node and k loops. The space H^ 1 ^ can be obtained from \ by adding to 



any graph from H^ k a new node t + 1 and k edges in the following k steps. At the ith step, we add one 
edge between the new node and one of the existing nodes 7. If 7 ^ t + 1, then it corresponds to a group of 
nodes 71, . . . , 7& in H^l +l ~ l \ The sum of degrees of 71, . . . , 7^ equals the degree of 7 in the graph before 
the ith step. We denote this degree by deg^. So 

deg. As) + k(a — 1) 

Prh — s) — ' - - 1 < s < t 

17 S) (a + l)(kt + i)-l' L ^ S ^ T - 

If 7 = t + 1, the corresponding group in H^f +l has only i — 1 nodes. Hence, 

deg tii (* + l) + (z-l)(a-l) + a 



Pr( 7 = t + l) 



(a + l)(kt + i) - 1 



We want to express any value r(d,t) in terms of some values with smaller t. Let us consider the 
transition from H^' k to H^ 1 ^. Let r(d, t, i) denote the average number of nodes of degree d, not including 
the last node t + 1, before the ith step, and r(d, t, i + 1) — the similar number after the ith step. Let 7 
be a head of the edge added in the ith step. Then, 

t t 
r(d,t,i+l) = J]Pr(deg M+1 ( S ) = d) = ^(Pr(deg M+1 (s) = d, 7 = s)+ 

s=l 8=1 

i 

+ Pr(deg M+1 (s) = d, 7 ^ s)) = ^(Pr(deg M (s) = d - 1, 7 = s)+ 

8=1 

+ Pr(deg tji ( S ) = d, 7 ^ ,)) = g fpr(deg tii ( S ) = d - 1) ( ( ^ ~ ^y^- + 



+ Pr(deg M (s) = d) 



d + fc(a - 1) 
(a + l)(H + i) - 1 



/, , .x (d - 1) + k(a - 1) ., / d + Ha-1) \ 

= rfd-l,t,i)7 rr, — ^ - + r(d,t,i) [l--, 777 f- — ■ (2) 

v ' ' J (a + l)(kt + i) - 1 v ' ' ' \ (a+l)(kt + i)-l) KJ 

By definition, 

r(d,t) = r(d,t,l), r(d, t + 1) = r(d, t, k + 1) + Pr(deg tijfc+1 (t + 1) = d). (3) 

The function r(d,t) is completely determined by the equations (j2J), © and the starting condition 

r{d,l) = {d = 2k}. (4) 

The equation (jHJ) includes the function Pr(deg t + 1) = d). Obviously, 

Pr(deg 4jfc+1 (t + 1) = d) = 0, d < fc or d > 2k. (5) 

The minimal value deg tfc+1 (t + 1) = k is obtained when no one of the edges is a loop. In this case, 
deg ti (t + 1) = i — 1 for all i, so 

pr ( de gl , i+l(t + 1) = t) = n (1 - (a+1)( :; +1) _ 1 ) - 1 + O (I) . 

(Note that a constant in 0() depends on a and fc). 

2k 

Because ^ Pr(deg tfc+1 = d) = 1 and Pr(deg tfc+1 = d) ^ 0, we get 



d=k 



Pr(deg tjfc+1 (t + 1) = d) = O , k < d ^2k. 



Since d is bounded in the last equality, its right hand side can be equivalently written as O (^r). 
Let 



c(d) 



B(d-k+ka,a+2) j ^ i 
B(ka,a+1) ' " ^ ft > 



0. 



d < k. 



(6) 



If d > fc, then 
c(d-l) B(d - 1 - k + ka,a + 2) 



c(d) 



B(d-k + ka,a + 2) 



r(d — 1 — + ka) /r(d + 1 — A; + &;a + a) d + 1 — &; + /ca + a 



Also 



c(fc) 



r(d - fc + ka)/T(d - k + ka + a + 2) 
B(fca,a + 2) r(a + 2)/T(ka + a + 2) a + 1 



d — 1 — + ka 



B(ka, a + 1) r(a + l)/T(/ca + a + 1) fca+a+1 

In particular, c(d — 1) > c(d), so c(d) < c(&) < 1 for all d ^ k. 

For the rest of the proof, we will assume that d ^ k. Note that of course this does not imply d— 1 ^ k. 
When d grows, the asymptotic behaviour of c(d) is 



In c(d) = In 



T(a + 2) T(d-k + ka) T(a + 2) 

B(fca, a + 1) r(d - fc + A;a + a + 2) ~~ B(ka, a + 1) 



+ (d - k + feo) (ln(d - k + fca) - 1) - (d - k + ka + a + 2) (ln(d - fc + A;a + a + 2) - 1) + O I - 

T(a + 2) . , , , w , , —k + ka „. ,, , , „ Wl , — A; + A;a + a + 2 
ln^rr+ — + (d- £; + fca)(lndH 1) - (d - k + fca + a + 2)(ln d H 1) 



B(£;a,a + 1) 



d 



d 



+ °'d 



c(d) 



T(a + 2) J _ 2 _ a 



B(£;a,a + 1) 



tT 



ln 



1 + 



T(a + 2) 
B(/ca, a + 1) 

1 

d 



(a + 2)lnd + - 
d 



(7) 



Let 



% I 

f(d, t, z) = r(d, t, i) — c(d) I t + - — 



k k(a + l) 

It is easy to see that the theorem is equivalent to f(d,t,i) = 0(1). Using (j2j), we obtain 



r (d — 1, £, z 



(d - 1) + Jfe(a - 1) 
(a + l)(kt + i) - 1 



+ r(d,t,i) 1 - 



d + fe(a — 1) 
(a + l)(H + z) - 1 



z 1 
r(d,t,i + l) - [ t + - - 



fc fc(a + l) 



c( d- 1)1^1+^4+ 



+ c(d) 1 - 



d + A;(a — 1) 
(a + l)(fct + i) - 1) 



r(d,t,z + l) - t + 



a+ l)(fct + i) - 1 
i 1 



fc fc(a + l) 



c(d)x 



(1 - [d = k])(d+ 1 - k + A:a + a) + (a + l)(kt + z) - 1 - (d + fc(a - 1)) 



r(d, i, z + 1) — c(d) 



(a + l)(fct + i) - 1 
[a + l)(kt + i + 1) - 1 - [d = k](l + ka + a) 
fc(a + 1) 



= f(d, t, z + 1) + [d = k]c(k) 1 + ka ~ { ~ a = r{d, t,i + l)+ ( 8 ) 
Let C — C(a, k) be a sufficiently large constant which will be determined later. We claim that 



f(d,t,i) + (z - 1) 



[d = k] 



k 



C 



min{l, ka} 



d + ka \ (a + 2)(t + l)(d + ka) 



i-1 



(9) 



for all i = 1, . . . , k + 1 and for all natural d ^ k and t. Note that this implies r(d,t,i) = O (^) and 
Theorem 1. 

The equations ()2]), ([3]), (j4]), fl5]) imply that r(d, t, i) — if d > kt + i — 1 + k. In this case, using ([7]), we 
obtain 

r(d, t, i) = -c(d) (t + i - ) = O = O (cT 1 -) , 

so if d > kt + i — 1 + k, © is true for all sufficiently large values of C. 

Now assume d^kt + i — l + k. We will prove ([9]) by induction on i and, for fixed t, on z. The basis of 
induction t = 1, . . . , 1 + |_^J and any i — 1, . . . , k + 1 obviously holds for all sufficiently large values of C. 

Now let t ^ 2 + L^J and let © hold for t - 1. Using ©, we obtain 

f (d, t, 1) = f (d, t - 1, k + 1) + Pr(deg t _ ljfc+1 (t) =d)= r(d, t - 1, k + 1) + [d = k) + O (J^J . 
Therefore, 

g / min{l,fco} \ fc of— ^ < 

|r( ' ' )! " T+ka~ V ~ {a + 2){t + l)(d + ka) ) + \cPi) " 

^ g Cmin{l,fco}/(o + 2) /J_ 
^ d + ka {t + l)(d + ka) 2 + \~cPt 

Thus, induction step on t is proved for all sufficiently large values of C . 

Finally, let t ^ 2 + |_^J > 1 + ^, z > 1 and let (El) hold for i — 1. We temporarily denote T 
(a + + z — 1) — 1. Note that T depends on i and z, but not on d. From (jSj) we obtain 

r W t, i) + (i - 1)1^1 = (i- 2)11^1+ 

+ f(rf - 1, t, i - +*(«-!) + ^ M _ 1} A i + *(. - 1) 



f(d-l,t,i-l) + (i-2) 



[d- 1 = k)\ (d-l) + k(a- 1) 



k J T 



+ (w, i - 1) + (i - 2)1^1) (l - d + + [fc < i < * + 1]0 (i) . 

The remainder term [A; ^ d ^ /c + 1]0 (|) can be written as O (^)- The assumptions d^kt + i — 1 + k 
and t > 1 + ^ imply that 1 - (a+ ^ t ( °7_l)-i ^ °- If d > At, then 

C (d - 1) + Jfe(a - 1) C A d + fc(a-l) 



d — 1 + ka T d + ka \ T 

C ( If d+ka 



d + ka \ T \ d — 1 + ka 
C ( x _ 1 f _ {d - 1) + k(a - 1)\ \ C 



((d - 1) + k(a -l))+d+ k(a - 1) 



d + ka \ T \ d — 1 + ka J J d + ka \ T(d — 1 + ka) 



^ C / A;min{l,A;a} 



If d = k, then 

C ^ d + A;(a-l)\ C ( ka(k + ka)\ < C / Armin{l,A;a} 



d + A;a\ T y d + A;a\ T(d+ka) J d + ka\ T(d + ka) 



In both cases, 



k 



^ C f min{l,A;a} V 2 / fcmin{l,fca} 

^ d + ka\ ~ (a + 2)(t+l)(d + ka) J \ T(d + ka) 



+ o[±- t ). (10) 



Note that 

fcmin{l,A;a} / k min{l, ka} 



1 



((a + l)(kt + i-l)-l)(d + ka) \ (a + 2)(kt + k)(d + ka) 



fcmin{l,fca} Q (\ 



(a + l)(a + 2)kt(d + ka) \ \t 

where the left part is strictly less than the first term in the right part, so that the second term in the right 
part is positive. Now ( ITUj) implies for all sufficiently large values of C. 
Theorem 1 is proved. 

4 Proof of Theorem 2 

By definition and linearity of expectation, 

cov(i*(di, t), R(d 2 , t)) = E(R(di, t)R(d 2 , t)) - r(d u t)r(d 2 , t) = 

t 

= E ^ [ de S s i = di,degs 2 = <k] - r(d 1 ,t)r(d 2 ,t) = 

S1,S2=1 

= yj Pr(degsi = di,degs 2 = d 2 ) + [d\ = d 2 )r(di,t) - r(di,t)r(d 2 ,t). (11) 

Slj=S2 

We will estimate the sum 



r 2 (d 1 ,d 2 ,t)= ^ Pr(degsi = d 1 ,degs 2 = <k) 

Sl,S 2 = l 

as we have done it for the function r(d, t) in the proof of Theorem 1. 

As with r(d, t), we define a function r 2 (di, d 2 , t, i) as the value of r 2 (di, d 2 , t) before the ith step in the 
transition from H^ k to H^ l \ The recurrent equation is deduced similarly to (T2]). For fixed s\ and s 2 , 
there are three non-intersecting cases: 7 = Si, 7 = s 2 , and 7 {si, s 2 }. In the first case, we get 

Pr ( d eg M+ i(si) = di,deg M+1 (s 2 ) = d 2 ,j = s x ) = 

= Pr(deg t>i (si) = di - l,deg M+1 (s 2 ) = 4,7 = «i) = 

di - 1 + fc(a - 1) 



Pr (degi,i(si) = di - l,deg f)i+1 (s 2 ) = d 2 ) 



(a + l)(kt + i) - 1' 

The second case is the same with d\ and <i 2 interchanged. In the third case, we get 

Pr (deg M+1 (s!) = d!,deg M+1 (s 2 ) = d 2 ,7 ^ «i,7 ^ «a) = 

= Pr(deg M (si) = di,deg M+1 (s 2 ) = d 2 ,7 ^ s 1; 7^ s 2 ) = 

, , s , , ( \ di + k{a-l) + d 2 + k{a-l) \ 
= Pr(deg tii ( Sl ) = dl , deg M+1 ( S2 ) = d 2 ) (fl + 1)(fct + _ 1 J > 



so the final formula is 

(di - 1) + k(a-l) 



r 2 (di,d 2 ,t,i + 1) = r 2 (di - l,d 2 ,t,i) 



(a + l)(kt + i) - 1 



, , ± .Ad 2 - 1) + k(a-l) 

+ r 2 d 1; d 2 - 1, t, « r-j4rrrT + 

(a + l)(kt + i) — 1 



+ ^.*.M)(l- ^Vl '- <12) 



By definition, 

r 2 (d 1 ,d 2 ,t) = r 2 (d 1 ,d 2 ,t, 1), 



r 2 (di, d 2 , £ + 1) = r 2 (di, d 2 , t, fc + 1)+ 

t 

+ Pr ( d eg t , fc +i(s) = di, deg tjfc+1 (t + 1) = d 2 ) + 

8=1 

t 

+ ^Pr(deg tjfc+1 (s) = d 2 ,deg tfc+1 (t + 1) = di), (13) 

8 = 1 

and the starting condition is 

r 2 (di,d 2 , 1) = 0. 

The equation (JT2j) includes a function 

t 

r 2 (d 1; d 2 , *, «) = Pr ( de §M( s ) = d h deg M (t) = d 2 ) (14) 

s=l 

(and the same function with swapped arguments), so we will first estimate r' 2 . Again, we write a recurrent 
equation. For fixed s, there are three non-intersecting cases: 7 = s, 7 = t+1, 7 ^ {s, t+1}. If deg t i+1 (s) = 
d%, then deg t j(s) equals di — 1 in the first case and di in the two other cases. If deg ti+1 (t + 1) = d 2 , then 
deg ti (i + 1) = d 2 — 2 in the second case and d 2 — 1 in the two other cases. Calculating probabilities, we 
obtain 

ifi ifi {di — 1) + k(a — 1) 

r 2 (di,d 2 ,M + 1) = Wi - Ma - — , ttttt , ^ — T+ 

(a + l)(fct + z) — 1 

, . , , , (d 2 - 2) + (i - l)(a- 1) + a 

+ r di, d 2 - 2, t, z l 2 r + 

(a + l)(fct + z) — 1 

/ di + fc(a - 1) + (rig - 1) + (i - l)(a - 1) + a \ 
+ r 2 (d 1; d 2 - 1, t, z) ^1 (a + l)(Jfct + i)-l J ' 

Before the 1st step, the node t has degree 0, so 

t 

r' 2 (d u d 2 ,t,l) = ^Pr(deg 4il (s) = d x )[d 2 = 0] = [d 2 = 0]r(d 1; t). 
s=i 

We continue to use notation ([6]). Obviously, r' 2 (di,d 2 ,t,i) = when d 2 > 2(z — 1) or d\ ^ 2(kt + i). If 
d x < 2(Art + i) and d 2 ^ 2(i-l), then ■ 0(d^ 2 t) = O^ " 1 ) = O(d^) and ^^-O(d^) = O(d^). 
Now it is easy to see that 

r 2 (d!,d 2 ,t,z) = [d 2 = i- l]c(di)t + 0(dr 1 ). (15) 

Since r 2 (di, d 2 , t, z) = when d 2 ^ 2/c, the remainder term in (1151) is zero when d 2 ^ 2A; and can be written 
as O 



did! 



3+1 



Since r(d\,t) = when d\ > kt + k, r' 2 (di, d 2 , t, i) = when d\ + d 2 > kt + + 2(i — 1). In particular, 
r^(d x , d 2 , t, A; + 1) + r 2 (rf 2 , di, t,k + l) = when d x + d 2 > kt + 2k. By ([12]) and (IT51) . r 2 (d x , d 2 , t, i) = 
when d\ + d 2 > kt + 2k + {i - 1). 

For the rest of the proof, we will assume that d\ ^ k and d 2 ^ k. Note that of course this does not 
imply di — 1 ^ /c and d 2 — 1 ^ k. 

Let 

f 2 (di, d 2 ,t,i) = r 2 (di, d 2 ,t,i) - c{di)c(d 2 ) ( t + \ — — -, : — J ( t + * 



k k(a + l)J \ k k(a+ 1) 

We temporarily denote T = kt + i — We express r 2 in terms of f 2 and use ( Fl2|) . In the expression, 
there are terms with r 2 with various arguments. Now we transform the terms without r 2 from the right 
part of PI). 

TT + l( ( dl -i) + fc ( -i) 

+ c(dx)c(<fc - l) ^" 1 ^^"^ + C (di)c(d,) f 1 dl + ^ + 2 * ( ° - 1} 



(a + l)T s A ' v ' y V (a + l)T 

I _c (dl) c(d 2 ) (a - [i = *]) (jrjTfjr + 

d 2 + l — k + ka + a d\ + d 2 + 2A;(a — 1) 



+ (1 - [d 2 = k\) 



(a + l)T (a+l)T 



_ TT + 1 / 2 (a + 1) + (a + 1)T - [di = fc](l + fca + a) - [d 2 = fc](l + A;a + a) 

" I^ c(dl)c(4) ^ (aTT)T 

T + 1T + 2 (T + l)c(d 2 ) (T + l) c (d0 
= "~ fc ^ c ( rf i) c ( rf 2) - [di = k] [d 2 = k] . 

The first term equals the term without r 2 in the left part, so 

~ i r t ■ \ . , .v (di — 1) + k(a — 1) 

r 2 (d 1 ,d 2 ,t,i + 1) = r 2 (di - l,d 2 ,t,2)- — -+ 

(a + + zj — 1 

, , .Jd 2 -l)+k(a-l) . , , . / di + G? 2 + 2k(a - 1) 
+ r 2 (d 1 ,d 2 -l,t,z V > nx \, ± + r 2 {d u d 2 ,t,t) 1 



(a + l)(Jfe* + i)-l (a+l)(fct + i) - 1 

(fct + z - 4t + !) c(d 2 ) (kt + i — T7 + l) c(di) 
- Mi = ifcp f 2 ] - [d 2 = k]- °g 

Relations (|T3|) and (|T5|) imply 

r 2 {di,d 2 ,t, 1) = f 2 (d u d 2 ,t- l,k + 1) + [d 2 = fc]c(di)t + [di = fc]c(d 2 )t + O ( +1 + ' 



Let 

r(d' 1 -k + ka) 



d^~^d2 did*^ J 
; d-^ ^ A^j d<2 ^ 



ci(d' 1; d 2 ) = { (d' 2 + k(a - l))T(d[ - k + ka + a + 1) 

0, d[ < k or d' 2 < fc 



By definition, for di > k, 



For d 2 > fc, 



Ci(di — 1, d 2 ) di — + ka + a 

Ci(di, d 2 ) d\ — k + ka — \ 

ci(di,d 2 — 1) d 2 + fc(a — 1) 

ci(di, d 2 ) d 2 — 1 + k(a — 1) 



Similar to (J7J), 

Moreover, 

ci(di,d 2 ) B(ka, a + 1) 



-l-a 



1 + 



d 2 + k(a — 1) 
T(di — k + ka)Y{d x -k + ka + a + 2) 



c{di) 



T(a + 2) (d 2 + k(a - l))T(di - k + ka + a + l)T(di - k + ka) 

B(ka, a + 1) di — k + ka + a + 1 



T(a + 2) d 2 + k{a - 1) 
Let C = C(a, k) be a sufficiently large constant which will be determined later. We claim that 



i — 1 

r 2 (d 1 ,d 2 ,t,i) H — ([d 2 = k]c(di) + [d\ = k]c(d 2 ))t 



< C(ci(di, d 2 ) + cx(d 2 , di)) kt + 



a+\ 
a + 1 



(16) 



for all i = 1, . . . , k + 1 and for all natural di ^ k,d 2 ^ k, t. Since both parts of ( TIB"]) are symmetric in di 
and d 2 , it is sufficient to consider the case d\ ^ d 2 . 

If d\ + d 2 > kt + 2k + (i — 1), then r 2 (di, d 2 , t, i) = and 



|f 2 (di, d 2 , t,i)\ _ c(di)c(d 2 ) (t + f - ^Tiy ) (t + - ^rry) 



ci(di, d 2 )t 



ci(di,rf 2 )* 



d 2 + k(a — 1) 
d\ — k + ka + a+1 



c(d 2 )t = O 



did\+ a J ' 



Since d 2 ^ ^i±^i ^ ^t+k, the right part is bounded. Obviously, [d 2 = k] = and [d\ = k]c(d 2 )t = O 

is bounded too. Thus ( [TBI holds when di + d 2 > kt + 2k + (i — 1) for all sufficiently large values of C. 

Let d\ + d 2 ^ kt + 2k + {i — 1). We will prove (JIB!) by induction on i and, for fixed £, on i. The basis of 
induction t = 1, . . . , 2 + |_^J and any i = 1, . . . , k + 1 obviously holds for all sufficiently large values of C . 

Let t ^ 3 + [_^J and let ( [TBI) hold for t — 1. We continue to use the restriction d\ ^ d 2 . Thus, 



|f 2 (di,d 2 ,t, 1)1 



f 2 (di, d 2 , £ - 1, + 1) + [d 2 = fc]c(di)t + [di = fc]c(d 2 )f + O 

a + | 



d? +1 d 2 



< C(ci(di, d 2 ) + ci(d 2 , di)) - 1) + 



Since c\{d\, d 2 ) = O ( ^+i d2 J > the right part is less than C{c\{dx, d 2 ) + ci(d 2 , di)) \ kt + ^+f-J for all suffi- 
ciently large values of C. This completes the induction on t. 

Let t > 3 + L^J > 2 + ^, i > 1 and let (JTSJ hold for i - 1. Then, 

i — 1 

f 2 (di, d 2 , t, z) H — ([d 2 = fc]c(di) + [di = k]c(d 2 ))t = 

. (di - 1) + k(a - 1) , . (d 2 - 1) + A;(a - 1) 

= f 2 (di - 1, d 2 , t, i - 1) . 1 . - + r 2 (di, d 2 -l,t,i- 1) , 1 _ w / . ^ \ + 



+r 2 (di,d 2 ,t,i - 1) ( 1 



a+ l)(fct + i - 1) - 1 

di + d 2 + 2A;(a - 1) 
~ (a + l)(fct + i - 1) - 1 



[di = jfe] 



a+ l)(ife* + z - 1) - 1 

(kt + i-1-^ + 1) c(d 2 ) 
k 2 



(kt + i - 1 - 4t + l) c(di ) z'-l 
■[d 2 = fc]- ^ ^-^ + — ([^2 = fc]c(di) + [di = k]c(d 2 ))t 



fr 2 (d 1 -l,d 2 ,t,i- 1) + — ^ ([d 2 = k]c(di-l) + [dt-1 = k]c(d 2 ))?j 

( i-2 
+ lr 2 {di,d 2 -l,t,i- 1) + — — ([d 2 -l = k]c{dx) + [d x = k)c(d 2 -l))t 

+ (f 2 (d u d 2 , t,i-l)+ ^^([d 2 = k)c(dx) + [di = k]c(d 2 ))t\ (l - - 



(di-1) + k(a - 1) 
(a + l)(H + z- 1) - 1 + 
(4-1) + fc(a- 1) 
(a + l)(jfet + i - 1) - 1 
di + d 2 + 2k(a - 1) 



+ 



[a + l)(kt + i - 1) - 1 
+ [d x < A; + l]0(c(d 2 )) + [d 2 ^ k + l]0(c(d 1 )). 

The assumptions d\ + d 2 ^ fct + 2 A; + (z — 1) and t > 2 + imply that 1 — {a+^^t+i-\)-\ ^ ^- ^ ^ ne 
induction hypothesis 



z — 1 

r 2 (d 1 ,d 2 ,t,i) H — ([d 2 = k]c(di) + [d\ = k]c(d 2 ))t 

K 



[kt + 



(z - 1) (ci(di - 1, d 2 ) + ci(d 2 , d x - 1)) 



o+j 
a + 1 

+ (ci(di, d 2 — 1) + ci(d 2 - 1, di)) 



(di - 1) + /c(o — 1) 



a + l)(fc£ + i - 1) - 1 
(d 2 - 1) + k(a - 1) 



+ (ci(di, d 2 ) + ci(d 2 , di)) 1 - 



a + l)(H + z - 1) - 1 
di + d 2 + 2£;(a - 1) 



(a + l)(kt + i- 1) - 1 

+ [dx ^ fc + l]0(c(d 2 )) + [d 2 ^ fc + l]0(c(di)). 



Since 



ci(di — 1, d 2 ) (di — 1) + k(a — 1) ci(di, d 2 — 1) (d 2 — 1) + fc(a — 1) 



ci(di,d 2 ) (o + l)(kt + i - 1) - 1 ci(di,d 2 ) (a + l)(fct + z - 1) - 1 

di + d 2 + 2Ha — 1) . r d\ — k + ka + a 

+ 1 - t — — , . 7 = (1 - [di = k\)- 



(a + l)(kt + i- 1) - 1 



(a+l)(kt + i-l)-l 



d 2 + k(a-l) d t + d 2 + 2k(a - 1) a 

+ (1 - [d 2 = k\)- 7VW~. . ■ — 7 + 1 - ~ TTTT. . ■ 7\ 7 ^ 1 + 



(a + l)(kt + i - 1) - 1 



{a + l)(kt + i - 1) - 1 



a + l)(kt + i- 1) - V 



we have 



z - 1 



r 2 (di,d 2 ,t,i) H — ([d 2 = fc]c(dj) + [di = fc]c(d 2 ))£ 



< C(ci(di, d 2 ) + ci(d 2 , di)) fci + 



a + 1 



•-1)1 + 



(a + l)(kt + i - 1) - 1 



+ 



+ [di < fc + l]0(c(d 2 )) + [d 2 ^ fc + l]O(c(d0). (17) 



Since 



kt + 



a + 



- \kt + 



a+1 

a + | 



a + 1 



fct + 



a + 1 



(a + l)(kt + z - 1) - I J a + 1 
1 



(a + l)(kt + i - 1) - 1 (a + l)(kt + i - 1) -1 



a + - ) (fct + i - 1)- 



Q + 2 

a + 1 



— afct — a 



a + 2 

a + 1 



(a + l)(Jfet + z - 1) - 1 V^" + a+1 



i-2) 



is always positive and tends to a nonzero constant limit as t grows, it is bounded from below by a positive 
constant. Therefore, for all sufficiently large values of C, the inequality ( IT7]) implies the inductive step by 
z, and so (fTSj) holds. 



As a consequence of (ITSl) . we obtain 

f 2 (rf 1 ,4,M) = o(^- + ^ TT ). 
The proven bound, the representation (TlTl) and Theorem 1 give the following bound: 

If di ^ d 2 , the maximum among the first three terms on the right-hand side is 0(<ij~ 2_a t); otherwise, the 
maximum is 0(d 2 2 ~ a t). This proves Theorem 2. 



5 Proof of Theorem 3 



We will use the notation N(si, s 2 ) for the number of edges between nodes si and s 2 . As usual, N tj i(si, s 2 ) 
is the value of N(s%, s 2 ) in the graph before the ith step. 
First, we define a function 



/ 



f(d 1 ,d 2 ,t,i) = E-, 



t t 



\ 



[degsi = di, degg 2 = d 2 ]Af(sx, s 2 ) 



V 



•1=1 S2=l 



/ 



It is easy to see that EX(di, d 2 , £) = /(di, d 2 , £, 1). 

Recurrent equations on / are deduced as it was done in the previous sections. The sum (ITSl) does not 
include the last node, so N(s±, s 2 ) does not change while adding a new edge. Thus, the ith step acts on / 
as in the case of r 2 (compare with (1121) ): 

/(d 1; d 2 , t, i + 1) = f(dt - 1, d 2 , t, i) {dl + 



(a + l)(fct + i) - 1 



,/ , , x ( rf 2 - 1) + Ha - 1) 
+ f(d 1 , d 2 - 1, t, i) \ 2 J„, \, ^+ 



(a + l)(H + i) - 1 

+ f(d h d 2 ,t,i) ( 1 - 



dx + d 2 + 2k(a - 1) 
(a + l)(kt + i) - 1 



(19) 



Second, we define a function 

d 2 , t, i) = E tji ([deg(t + 1) = d 2 ) ^[deg s = tfc] JV(f + 1, s) 

\ 8=1 

Obviously, 

/(c?i, d 2 , £+1,1) = /(di, d 2 , t, + 1) + p(c?i, d 2 , t,k+l) + g(d 2 , di, t,k+l) 
and since N(t + 1, s) =0 before adding any edges from the node t + 1, 

g(d 1 ,d 2 ,t, 1) = 0. 



(20) 



(21) 



(22) 



We now consider one summand of the sum ( 1201) and the ith step. Let the new edge link nodes t + 1 
and 7. We have three non- intersecting cases: 7 = s, 7 = t + 1, 7 ^ {s, t + 1}. Note that 

[7 = s, deg M+1 (t + 1) = d 2 , deg M+1 s = di]N t>i+1 (s) = 

= [7 = s,deg M (t + 1) = tfc- l,deg M s = d x - l){N tyi {s) + 1), 



[7 = t + 1, deg t i+1 (t + 1) = d 2 , deg ti+1 s = d x ]JV tji+1 (s) = 

= [ 7 = t + 1, deg M (t + 1) = d 2 - 2, deg M s = di]N t ,i(s), 



[7 £{s,t+ 1}, deg tji+1 (t + 1) = d 2 , deg M+1 s = £Z x ]iV tji+1 (s) = 

= [7 £ {s, t + 1}, deg tji (t + 1) = da - 1, deg M s = d x ]iV M (s). 

Taking the expectation and using the definition (TH1) . we obtain 

. (d x — 1) + k(a — 1) 

5f(di, d 2 , t, i + 1) = (</(di - 1, d 2 - 1, t, i J + r 2 {di - 1, d 2 - 1, t, «)b — , , — 7+ 

(a + l)[kt + 1) — 1 

/_, , o , .x(d2-2) + (i-l)(o-l) + a 

+ <? d x , d 2 - 2, t, z p — r — + 

(a + !)(«;£ + 2) — 1 

+ ^,4-1,^)^1 (a + l)(fct + 0-l )' (3) 

Third, we derive a bound on (7. Obviously, g(di, d 2 , i, z) = when d 2 > 2(z — 1) or di ^ 2(/ct + z). If 
d x < 2(kt + i) and d 2 <: 2(i - 1), then ■ 0(dr a_1 ) = O^f 1 ) = 0(r x ) and ^ • O^ 1 ) = Oft" 1 ). 
Remember that we have proved the bound (I15p on r' 2 . It is easy to see now that 

g(d u d 2 , t,i+l)=i[ik = z]c(d x - l) {d] [ ~ "^ffi" ~ l) + O (jj . (24) 



Finally, we are ready to study /. For the rest of the proof, we will assume that d x ^ k and d 2 ^ k 

T(a+2) / 1 -I \ r(fcq+aH 

B(fco,o+l) V a "T J-J r(fca) 



We denote A = B ^ a+ ^ , = (a + 1) wt^" 1 " , -Di = di — + ka, D 2 = d 2 — k + ka for brevity. By 



definition, 

c(d x ) = A r(A 



r(Dj + a + 2)' 
Let cx(di,d 2 ) be defined recurrently as follows: 

c x (k,k) = 0, 

(A-lXcxtrfx-l^ + cCd!-!)) 

c x {di,k) = — — — ■ — , di>k, 

Di + ka + a + 1 

, x P 2 -I)(c x (fc,d 2 -1)+C(d 2 -1)) 

C * (M2) = D 2 + ka + a + l ' d2>k > 

c x (d h d 2 ) = (A-l)cx(di-l,d a ) + (A,-l)cx(di,d a -D 4)d2> ^ 

D x + D 2 + a + 1 



Let 



c 2 (di,d 2 ) 



c 3 (d x ,d 2 ) 



r(Dj + 2)r(^ 2 + 2)rpi + d 2 + a + 2) ' 
r(D 1 )r( J D 2 )r( J D 1 + J D 2 + i) 



r(A + 1)T(D 2 + l)r(Di + D 2 + a + 2) ' 
Obviously, these functions are symmetric. If d x > k, 

c 2 (d!-l,d 2 ) = {D 1 + l)(D 1 + D 2 + a + l) 
c 2 (d x ,d 2 ) ~ (D 1 -l)(D 1 + D 2 + 2) ' 

c 3 (di - l,d 2 ) _ + £> 2 + a + 1) 

ca^ds) ~ (£>! - l)(Di + D 2 ) ' 

Thus, for di, d 2 > fc, 

(£>! - l) C2 (d! - 1, d 2 ) + (£> 2 - l)c 2 (d!, d 2 - 1) 



c 2 (di,d 2 ) 



D x + D 2 + a + 1 



( , A s {D l -l)c 3 {d 1 -l,d 2 ) + {D 2 -l)c 3 {d 1 ,d 2 -l) 

c 3 (d 1 ,d 2 ) = — 

L>i + L> 2 + a + 1 

c(d x ) r(D 1 + 2)r(to + 2)r(D 1 + fea + a+2) 



Afcac2(di, fc) fcar(Di + a + 2)r(fca)r(Di + fca + 3) 

A;a+1 r(L>i + 2)r(D 1 + A;a + a + 2) 
~~ D 1 + ka + 2 T(D 1 + a + 2)f(A + ka + 2)' 

Let (a) n = a(a + 1) . . . (a + n — 1) be the Pochhammer symbol. Let 2 F\ (a, 0; 7; z) = Xl^Lo ^"^j^" 
for —1, —2, ... be the hypergeometric function. According to pfl 15.1.1], if 7 — a — (3 > and 

I z I ^ 1, this series converges absolutely. We quote the following formula from pQ 15.1.20]: 

F^/wn r( 7 )r(7-a-/3) 

2-Pi(a, P, 7, 1) = ^7 rp7 

r(7 - a)r(7 - p) 

Thus 

r(D 1 + 2)r(D 1 + fc a + a + 2) ^ (a) n (fca)n 

rm.^x,)rma^x9) = sFl ' A. + fca + a + 2; 1) = 2^ 



r(Di + a + 2)V{D 1 + A;a + 2) iV ' ' ' 7 ^ (D x + A;a + a + 2) n n! ' 



Since all terms of the last series are positive and the first term is 1, 

c(di) ka + 1 



Akac 2 (di,k) Di + ka + 2 
Moreover, 

rpi + 2^^ + ka + a + 2) ^ (a) n+ i(A;a) n+ i 



rpi + a + 2)r(Di + ka + 2) ^ (A + to + a + 2) n+1 (n + 1 



_ 1 Q 2 fe (a + l) n (to + l) n < 

D 1 + ka + a + 2^(D 1 + ka + a + 3)Jn + 1)! ^ 

n=0 

a 2 A; 

^ 1 + — ; 2 Fi(a + l,ka+ 1; 2to + a + 3; 1) < 

L>i + to + a + 2 

< a 2 £; r(to + l)r(2to + a + 3) 

^ + D 1 + ka + l T(2ka + 2)T(ka + a + 2) ' 

c(di) < to + 1 / + to£ \ ^ _ T{ka + l)r(2to + a + 3) 



Akac 2 (d x , k) D 1 + ka + 2 \ A + + 1 / ' r(2to + 2)r(to + a + 2) ' 

In Theorem 3, we have three assertions. The first one says that EX(di,d 2 ,t) = cx{di,d 2 )t + 0^(1). 
The second one gives a bound for cx- The third one gives an asymptotic formula for cx- Now we shall 
show that our function cx admits the bound from the second assertion. This bound is equivalent to 

Aka ( c 2 (di, d 2 ) — ( 4 — — ) c 3 (di, d 2 ) J < c x (di, d 2 ) ^ Aka(c 2 (di, d 2 ) + Bc 3 (dx, d 2 )). (25) 



1 + ka 

To prove (|25|) . we use induction on di + d 2 . If d x = d 2 = k, the right-hand side of the inequality is obvious, 
and its left-hand side follows from 

c 2 (k,k) (2to + 2)(2to + l) Aka + 2 j 2 



c 3 (k,k) (ka + l)(ka + 1) ka + l ka + l 



\{ di > k and c?2 > & ; all the parts of f[2"5"j) satisfy the same recurrent equation, so ( I25p follows from the 
induction hypothesis. Due to symmetry, it remains to prove (|25p for d 2 — k, d\ > k. We have 



fca (c 2 (di,fc) - ^gc 3 (di,fc)) _ _ (A + l)(2 + 4fca) _ 
kac 2 (d 1 ,k) ~~~ (D 1 + ka + 2)(D 1 + ka + 1) " 

pi -fca) pi -fca + 1) 



pi + fca + l)pi + fca + 2)' 

In particular, fca (02(6?!, fc) — 2 X "|7^° 03(^1, A;) J > for d\ > k. Then, 

c x {di,k) 

Afca(c 2 (di,fc)-^c 3 (di,fc)) ~ 

pi + fca + l)pi + l)pi + ka + a + 1) c x (d u k) 



pi - fca)pi - fca + l)pi - 1) Akac 2 {d 1 - 1, k) 
pi + ka+ l)p! + 1) c x (di - 1, A;) + c(di - 1) 



pi - fca) pi - fca + 1) Afcac 2 (di - 1, k) 



> pi + fca + l)pi + l) /pi-1- fca) pi - fca) ^ fca + 1 



pi - fca) pi - fca + 1) V (A + fca) pi + fca + 1) Di + ka + 1 



Di + 1 {D 1 -l-ka ka + 1 



£>i — fca + 1 \ Di + fca D\ — ka 
Furthermore, 

D x - 1 - fca fca + 1 £>i - fca + 1 fca + 1 (fca + 1) A - fc 2 a 2 + 2fca + 1 
D 1 + ka Di - fca Z?i + 1 " D x -ka pi + l)p x + fca) ^ 

/ 1 A + fca + 1 \ _ , L>i + fca + 2 

^ (fca + 1} V^fc^ " p 1 + l)p 1 + fca) j=^ + 1)fca p 1 + l)p 1 + fca)p 1 -fca) > °' 

This proves the left hand-side of the inequality ( 1251) . 

Now, 

cx{di,k) 



Aka(c 2 (di, k) + Bc 3 (di, fc)) 



, , D pi + l)(fca + l) \ 1 pi+ l)pi + fca + a + l) c x {d u k) 
1 + ^77; ; tttt; ; ttt — tt: tttt; ; ^ 77 n : — rr ^ 



pi + fca + l)pi + fca + 2)y pi-l)pi + fca + 2) Afcac 2 (di - 1, fc) 

pi + fca + l)pi + l) 

^ pi + fca + l)pi + fca + 2) + B(D 1 + l)(fca + 1) X 

Di(fca + 1) fca + 1 A fcaS 



X 1 + ^ pi + fca) pi + fca + 1) ' L>i + fca + 1 \^ ' Dx + ka 



Furthermore, 



DAka + l) fca + 1 / kaB 
1 + B— -++- — + 1 + 



pi + fca) pi + fca + 1) D x + fca + 1 \ Di + ka 



D L + fca + 2 fca+1 \ / fca + 1 fca + 1 ' , 

+ B— -——) = [— — — ■ - — ■ — I < 0. 



Di + 1 Di + fca+iy \Dx + ka + l £>i + 1 

This proves the right-hand side of the inequality ( 125*1) . 

The third assertion of Theorem 3 (i.e., the asymptotic formula for cx{di,d 2 )) is derived from (1251) in 
the same way as a similar formula was derived from (j7j). 

It remains to prove the first assertion. We will use the bound cx(d\, d 2 ) = O ( ^-^p 

Let 

/ i 1 
f{di,d 2 , t, i) = f{di, d 2 , t, i) - c x {di, d 2 ) [t + 



fc fc(a+ 1)/ ' 



Let C = C(a, k) be a sufficiently large constant which will be determined later. We claim that 



(r , fl (D 2 -l)c(d 2 -l) r , ,i(-Di-l)c(di-l) 
f(d 1 ,d 2 ,t,i) + (i- 1) [rf 1 = A;] v \; ^ + [d 2 = fc]^ ] 



(a + l)fc 



(a + 1)& 
€ C | 1 - 



(a + 2)fc(f + 1) 



i-l 



• (26) 



Since r' 2 {d\, d 2 , t, i) = when di + d 2 > kt + k + 2(i — 1), (|22|) and (|23|) imply that g(di,d2,t,i) = 
when G?! + d 2 > H + + 2(i — 1). Consequently, (|T9|) and (I2T]) imply that f(di,d 2 ,t,i) = when 
d x + d 2 > kt + 2k + {i- 1). 

If di + cf 2 > fct + 2A; + (i - 1), then 



f(d 1 ,d 2 ,t,i) = -c x (d 1 ,d 2 ) lt + - 



k k(a + l) 



O 



(di + d 2 f- a 



so f[2"ET) holds for all sufficiently large values of C. 

Now assume di + d 2 ^ fci + 2/c + (i — 1). We will prove ( 1261) by induction on i and, for fixed t, on i. 
The basis of induction t = 1, . . . , 2 + |_^J and any i = 1, . . . , k + 1 obviously holds for all sufficiently large 
values of C. 

Now let t ^ 3 + L^J an d let (|26|) hold for t — 1. Since (1261) is trivial for d 1 = d 2 = k and symmetrical, 
we may assume that d\ > k. From (12 ip and (El]), we obtain 



f(d 1 ,d 2 ,t, 1) 



/(di, d 2 , i - 1, A: + 1) + k[d 2 = k\c{d x - 1) t^t^ + O (j 



^ C 1 



(a + 2)fci 



The right-hand side is less than C for all sufficiently large values of C. 

Finally, let t ^ 3 + [^J > 2 + i, i > 1 and let O hold for i-l. We reuse the notation T 
(a + l)(/c£ + i — 1) — 1 and again assume d\ > k. From (1191) we obtain 

/(di -l,d 2 ,t,i- 1)^-^ + [d 2 > fe]/(di, d 2 - 1, t, i - 1)^-^+ 



T 



T 



f f(d 1 ,d 2 ,t,i — 1) ( 1 — J ° 1 1, ^ 2 ) = f(d 1 ,d 2 ,t,i) + c x (d 1 ,d 2 ) ( - — ^-ttt + ~ I 



c x (di - 1, d 2 )(L>! - 1) + [d 2 > fc]c x (di, d 2 - 1){D 2 - 1) + c x (di, d 2 )(T - (£>j + D 2 )) 



(a + l)fc 



/(di,d 2 ,t,z) + [d 2 = fc]c(d 1 -l) 



(a + 



The assumptions di + d 2 ^ fct + 2A; + (i - 1) and t > 2 + i imply that 1 - = 1 - ^(hS-im ^ °' 

Since(D 1 -2)c(di-2) = [d a >Jfe+l]pi + a)c(di-l), 

(, _ 2)[d2 = ^-2)44-2)^-1 + (j _ 2)[d2 = t] (C 1T lW-l) ( t _ A + C 2 



(a + l)Jfe 



(a + l)fc 



T 



, . (L>i - l)c(di - 1) /. , , nJ Di + a D l + D 2 

(i-2)[d 2 = fc] 1 . ; M [di > A; + l]^^ + l 



(a + l)fc 



T 



Ifd 2 > fc,then^± + ^± + l-^±^ = 1-f < 1-i Ifd 2 = k, then^i + 1-^^ = ^ 
Thus, 

(Di - l)c(d! - 1) 



f(di,d 2 ,t,i) + (i - l)[d 2 = fe]- 



(a + l)Jfe 



For all sufficiently large values of C, the right-hand side is less than C ^1 — ( a+2 yk(t+i) ) ' so ^ ne induction 
on i is complete. 

Theorem 3 follows from the proven bound ( 1261) . 



6 Proof of Theorem 4 

We use the Azuma-Hoeffding inequality. 

Theorem 5. /73] / Let (X s )™ =0 be a martingale with \X s+ i — X s \ ^5 for s — 0, . . . , n — 1, and x > 0. 
Then 

P(\X n -X \ ^x)^2exp 

We fix di,d 2} t and denote X = X(di,d 2 ,t). Let G be a random graph in H^' k ; it has kt edges, sorted 
by the creation time. Let be a graph with s first edges. Let X s = E(X\G^), s = 0, . . . , kt. In this 
sequence Xo = EX, X^ t = X. By definition of the probabilistic space, the sequence X s is a martingale. 
We will estimate possible differences between adjacent elements of the sequence. 

We fix any s from to kt — 1. Let v be the head of the last edge in G^ s+1 \ so v is a random quantity 
depending on G. By definition 

X s = £ 7 Pr(t; = 7) E(X\G ( - S \v = 1 ), 
X s+1 = E(X\G^\v = v(G^)), 

where the sum is over all nodes of G. Hence it is clear that 

mmE(X\G (s) ,v = 7) ^ X S ,X S+1 ^ max E(X\G (s) , v = 7), 

7 7 

\X S -X S+1 \ ^ maxE(X\G {s \v = 7) - mmE(X\G {s \ v = 7). 

7 7 

Let 71 G argmin E{X\G^ s \v = 7) and 72 G argmax = 7). It is sufficient to prove an upper 

bound for 

E{X\G {s \v = 72) - E(X\G<-'\v = 71). 

We consider the sum 

t t 

X=Y.Y1 t de S Sl = rf i' de S s 2 = d 2 ]N( Sl ,s 2 ). (27) 

Sl=l «2=1 

Replacing the condition t> = 71 by the condition v = j 2 changes distributions of degrees of 7« and 
distributions of Nfy, *) = N(*,ji)] distributions of other values of N do not change. Thus distributions 
of all terms in the sum (127]) except those with {71, 72} H {si, s 2 } 7^ are the same for t> = 71 and v = j 2 . 
Let 

t t 

X' = ^2 S [degsi = rfi,degs 2 = c/ 2 ]X(si,s 2 ). 

81=1 S2=l 

{«i,aa}n{7i,7a}^0 

Then 

E(X - X'\G {s \ v = 71) = £(X - w = 72). 




Obviously, X' ^ 0. We have 

t t 
X' < J^[degsi = tZi,deg7i = d 2 ]N(s 1: 71) + J^fdegsi = di,deg7 2 = d 2 ]-W(si, 72) + 

Sl=l «l = l 

t t 
+ J^[deg7i = di,degs 2 = d 2 \N{^ u s 2 ) + J^[deg7 2 = di,degs 2 = d 2 }N(>y 2 , s 2 ) < 

S2 = l S2=l 

t t 

< [deg7x = d 2 \ ^ N ( s i,li) + [deg7 2 = d 2 \ ^ N(si,7 2 )+ 

Sl = l Sl = l 

t t 
+ [deg7i = di] ^2 N (7i, s 2 ) + [deg7 2 = di] #(72, s 2 ) = 

S2 = l S2 = l 

= [deg7i = d 2 ]d 2 + [degj 2 = d 2 ]d 2 + [deg7i = d^di + [deg7 2 = d^dx ^ 2(di + d 2 ). 

^ u = 7l ),S(X'|GW,u = 72) < 2(dx + d 2 ), 

\E(X'\G {s \v = 71) - £(X'|G< a \ u = 72)! < 2(d 1 + d 2 ), 
\X S - X s+1 \ < 2(di + d 2 ). 

Consequently, the sequence (X s ) satisfies the condition of Theorem 5 with n = kt and 5 = 2(di + d 2 ). 
Substituting x = c(di + d 2 )Vki in Theorem 5, we obtain Theorem 4. 



Thus, 
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